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This  report  presents  the  results  of  a preliminary  investigation  on 
the  stability  of  a nominally-two-dlmensional  laminar  boundary  layer  flow 
approaching  the  separation  point  with  respect  to  3-D  periodic  spanwise 
disturbances,  under  the  hypothesis  that  a steady  state  is  reached. 

The  basic  equations  of  the  disturbances  are  examined  with  an  assumed 
form  of  the  3-D  perturbations,  resulting  in  a system  of  ordinary  differen- 
tial equations.  Together  with  the  boundary  conditions  requiring  the  dis- 
turbances to  be  zero  at  the  wall  and  to  vanish,  asymptotically  at  infinity, 
a two-point  eigenvalue  problem  was  formulated.  Meksyn's  method  was 
adapted  for  the  calculation  of  the  basic  two-dimensional  flow  parameters 
entering  the  equations  in  the  region  near  the  separation  point.  This 
method,  which  is  superior  for  practical  purposes  to  Goldstein's  well-known 
analysis  of  separation,  requires  the  experimentally  determined  streamwise 
variation  of  the  pressure  at  the  edge  of  the  boundary  layer. 
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NOMENCLATURE 


also,  speed  of  sound 


D = 25(3/H)  - n(l  + AQ  + TQ)(3/3n) 

a.  b » c . coefficients  defined  in  Chapter  3 and  Appendixes  B and  C 

n n n n 
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e.n 


transformed  stream  function;  * J — dn 

o ue 

total  enthalpy  ratio,  HQ/He 

total  and  static  enthalpy  (H  = h + u2/2) 
streamwise  characteristic  length 
static  pressure 

coefficients  defined  in  Chapter  3 and  Appendix  B 
absolute  static  temperature 

streamwise,  normal,  and  sidewash  velocity  components 
streamwise,  normal,  and  lateral  coordinate,  respectively 
wave  number  (2iT/x) 
specific  heat  ratio , Cp/c.* 
boundary  layer  thickness 
transformed  Levy-Lees  coordinates 
d In  u_ 


d In  c 


= A 


o h. 


spanwise  disturbance  wavelength 
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y laminar  coefficient  of  viscosity 

v kinematic  viscosity  (u/p) 

p density 

T shear  stress; 

also,  t = / fQU»n)  dn 
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effective  edge  of  the  boundary  layer 
separation  point 
wall  surface 
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As  established  by  Gortler*,  a two-dimensional  flow  on  a concave  wall 
can  develop,  under  certain  conditions,  a system  of  3-D  time-dependent 
periodic  spanwise  disturbances  featuring  a row  of  alternating  streamwise 
vortex  pairs  superimposed  on  the  basic  two-dimensional  boundary-layer 
(see  Figure  1 ).  The  concavity  of  the  streamlines  and  the 

associated  normal  pressure  variation  are  the  primary  mechanisms  causing 
this  type  of  disturbance  flow. 

A similar  type  of  situation  has  been  found  to  occur  in  the  case  of  a 

two-dimensional  flow  approaching  a stagnation  point,  which  also  possesses 

o 

strongly  concave  streamlines  (Fig.  2).  In  addition  to  Gdrtler's  approach 
a series  of  papers  4*  5’  7 have  demonstrated  that  this  type 

of  disturbance  flow  can  exist  even  in  a steady  (i.e.,  neutrally  stable) 
state.  A variety  of  experimental  evidence  in  support  of  this  is  presented 
in  the  cited  papers. 

The  facts  presented  above  naturally  lead  to  the  idea  that  these  dis- 
turbance vortices  may  in  fact  occur  whenever  the  basic  two-dimensional 
flow  exhibits  sufficently  concave  streamlines.  In  particular,  this  might 
be  the  case  for  two-dimensional  separating  boundary  layer  flows  owing  to 
the  significant  streamline  curvature  that  develops  approaching  separation. 
In  fact,  a recent  experimental"  investigation  by  Settles,  et  al?,  featuring 


a detailed  examination  of  the  surface  flow  pattern  at  separation  by  the 
kerosene-graphite  technique  (see  Fig.  3)  exhibits  a distinctly  wavy  separa- 
tion line  that  suggests  the  presence  of  a spanwise-periodic  disturbance 
field.  Stimulated  by  these  observations,  the  present  authors  therefore 
undertook  a theoretical  study  of  possible  vortex-like  steady  state  distur- 


bances near  separation. 


We  consider  that  the  same  major  qualitative  properties  of  3-dimensional 


disturbances  in  both  the  flow  along  a concave  wall  and  in  a stagnation 
region  are  also  valid  in  the  case  of  a separation  region.  Consequently,  we 
may  assume  that  the  three-dimentional  disturbances 

(a)  are  periodic  spanwise  with  a wavelength  X that  is  larger  that  the 
boundary  layer  thickness  (accordi,  to  Gortler  X « 2-46), 

(b)  are  concentrated  in  a thin  layer  near  the  wall  whose  thickness  6' 
is  not  necessarily  equal  to  but  of  the  same  order  as  that  of  the  basic  bound- 
ary layer, 

(c)  have  a vortex  structure  featuring  a row  of  streamwise  vortex-pairs, 
and 

(d)  have  a normal  variation  of  the  pressure  which  is  a significant 
aspect  of  the  disturbance  mechanism. 

In  addition,  the  fact  that  the  separation  line  is  wavy  implies  that  the 
streamwise  disturbance  velocity  u‘  cannot  be  considered  zero  near  the  two- 
dimensional  separation  point.  Unfortunately,  the  experimental  data  available 
until  now  do  not  supply  any  indication  about  the  streamwise  variation  of 
disturbances. 

However,  taking  into  account  the  fact  that  ine  concavity  of  the  stream- 
lines is  likely  the  main  cause  of  this  type  of  flow  and  that  the  concavity 
has  a maximum  near  the  separation  line,  we  can  consider  that  the  streamwise 
variation  of  disturbances  is  small  enough  to  be  neglected  up  to  the  second 
order. 

The  model  of  the  flow  corresponding  to  the  above  mentioned  properties 
is  similar  with  that  of  Gortler,  being  represented  in  Fig.  4. 

This  report  describes  the  preliminary  results  of  the  theoretical  investi- 


gations on  the  spanwise  disturbed  2-D  flow  near  separation  location.  In 


a theoretical  model  followed  by  a detailed  formulation  of  the  appropriate 
linearized  equations  and  eigenvalue  problem  governing  the  small  distur- 
bance field.  Section  3 discusses  at  some  length  the  analytical  treatment 
of  the  basic  2-D  separating  boundary  layer  flow  itself;  we  show  that  this 
important  aspect  of  the  problem  is  best  treated  in  practice  by  an  adaption 
of  Meksyn’s  work  rather  than  using  Goldstein's  well-known  mathematical 
theory  of  separation.  Finally,  in  Section  4,  we  outline  the  suggested 
future  approach  to  the  solution  of  the  equations  developed  herein. 


2.  MATHEMATICAL  FORMULATION 


2.1  Preliminary  Appraisal  of  the  Disturbance  Mechanism 

Prior  to  a detailed  analysis  of  the  disturbance  field  equations,  it 
is  well  to  first  establish  the  physical  plausibility  of  the  hypothesized 
streamwise  vortex  disturbance  mechanism.  This  can  be  done  in  a very 
approximate  way  by  a local  application  of  Gortler's  instability  theory 
for  concave  streamline  flows;  the  proposed  mechanism  is  likely  to  occur 
if  under  typical  conditions  it  appreciably  exceeds  his  instability 
criterion. 

According  to  Gortler!  a concave  boundary  layer  flow  having  an  aver- 
agelongitudinal  radius  of  curvature  R will  be  susceptible  to  amplified 
streamwise  vortex  disturbances  whenever  the  following  instability  con- 
dition is  met: 

<eVL,V‘(j|L)  * .*  o, 

where  0*  is  the  momentum  thickness  and  L is  the  streanwise  distance  to 
the  local  curvature  region  under  consideration.  Now  to  apply  Eq.  (1)  to 
the  present  problem+,  we  must  supply  an  estimate  of  R typical  of  the  boun- 
dary layer  flow  near  the  surface  approaching  separation.  This  can  be 
found  from  the  following  analytical  expression  (derived  in  Appendix  A-l) 
that  relates  R to  the  observable  physical  properties  of  the  flow  (see  Fig.  4): 

R s:  (3p/^x)$/(  S2tw/^x2)s  (2) 

+ Since  the  boundary  layer  profile  upon  which  Eq . (1 )is  based  has  a signi- 
ficant wall  shear,  its  application  near  separation  is  only  a rough  approxi- 
mation because  it  is  not  known  how  the  RHS  is  influenced  by  an  adverse 
pressure  gradient. 
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We  thus  require  data  on  the  streamwise  pressure  gradient  and  second 


derivative  of  the  streamwise  wall  shear  distribution  [i.e.,  the  curva- 
ture of  the  graph  tw( x) ] at  the  separation  point. 

Fortunately,  there  exists  some  relatively  detailed  experimental  data 
at  conditions  close  to  those  pertaining  to  Fig.  3 from  which  these  quanti- 
ties can  be  readily  extracted:  Figures  5 & 6 illustrate  the  pressure  and 

wall  shear  distributions  measured  by  Sfei r for  a M = 2.64,  Re.  = 1.4  x 

CO  L 

5 

10  laminar  flow  separating  in  front  of  a small  angle  ramp  compression 
corner.  From  these  curves,  we  obtain  the  estimates 

■ .4  ,u>ref 

and 

ldtp/PJ/d(X/L)]s  - 1.7 

where  tw  ^ = .332  p e^e^/j^x  is  the  Blasius  reference  value  at 

X/L  - .25  with  L = Xj  = 2.7  inches  and  es*  = .07  inches.  Consequently 
we  find  by  substitution  into  Eq.(2)  that 

f “ 2.4^^L  = 129  (3) 

L ? 

M c 

oo 

which  in  turn  from  Eq.(l)yields  the  stability  parameter  value 

(e*UJv)^e*/R  » 28.6.  (4) 

Since  this  exceeds  the  critical  value  by  two  orders  of  magnitude,  it  is 
likely  that  the  streamline  curvature  approaching  separation  in  these 
typical  conditions  is  more  than  sufficient  to  give  birth  to  streamwise 
vortex-like  disturbances. 

2.2  Basic  Assumptions,  Order  of  Magnitude  Analysis  and  ueneral 
Disturbance  Equations 

In  order  to  derive  the  basic  disturbance  equations,  we  start  from  the 


5 


complete  Navier  Stokes  equations  of  viscous  compressible  flow  and  as  a 
first  step,  make  use  of  the  following  assumptions. 


I 


1.  The  flow  is  composed  of  a basic  two-dimensional  flow  uQ,  vq,  pQ, 
p , To  plus  a system  of  three-dimensional  steady  state  perturbations  u', 
v' , w* , p* , p' , T' . We  note  in  this  connection  that  the  basic  flow  is  not 
assumed  to  be  a parallel  shear  flow  (indeed,  in  the  present  study  such  an 
assumption  would  eliminate  the  disturbances  of  interest). 

2.  The  perturbations  are  small,  so  that  we  can  neglect  second  and 
higher  order  terms  as  compared  with  those  of  first  order  (linearization). 

3.  Provided  that  the  Reynolds  number  is  sufficiently  large,  the 
boundary  layer-like  approximations  apply  to  the  basic  flow  parameters, 
these  furnishing  the  well-known  estimations  that 

!o  _1_  _ / . . , 

*o  " Vie  IfPo^o 


V - -£u  * -r-  U , 

o £0  0 lo  " 


a « a 

5^  K<Sy- 

4.  Under  the  same  condition,  the  viscous  effects  on  the  disturbances 
are  concentrated  in  a thin  layer  whose  thickness  6'  is  considered  of  the 
same  order  as  6Q  and  small  as  compared  with  both  the  spanwise  and  stream- 
wise  characteristic  lengths.  Taking  into  account  that  Gortler's  value  for 
the  spanwise  wavelength  is  (2  - 4)  6Q,  we  assume  6'  /A  <<  1 and  thus 


a? 


« £ 


V 


As  for  the  streamwise  characteristic  length,  we  will  consider  S'/t'  <<  1 

and  hence  d/dx  <<  d/dy  (at  the  same  time  we  must  also  have  \/i'&  <5 '/A, 
A2/*.'2  « 1). 
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5.  As  in  Gortler's  flow  along  concave  walls,  the  streamwise  dis- 
ince  vorticity  w'x  dc 
This  can  be  expressed  by 


turbance  vorticity  w‘  dominates  the  normal  disturbance  vorticity  u*  . 

x y 


au' 

aw1 

aw1  a v 

dz 

■ ar  I 

< 

ar  ' "ar 

(5) 

6.  Considering  that  the  compressibility  does  not  change  essentially 
the  proportions  between  u',  v',  and  w',  we  obtain  the  following  estimation 
by  using  the  incompressible  disturbance  continuity  equation: 

v'  < 6' 


aw 

i f i 

au1 

ar 

+ 6 

ar 

(6) 


On  the  other  hand. 
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dzdx 


«: 


du' 


dz 


5 1 


a2W 


aZ‘ 


« 
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Consequently,  Eq.  (5)  gives: 


and  hence  from  Eq.  (6) 
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(7) 


(8) 


Applying  the  above  assumptions  and  assuming  also  a unit  Prandtl  num- 
ber and  Chapmann's  law  for  the  laminar  viscosity  (u/T  = const.),  the 
following  general  disturbance  equations  are  obtained: 

Continui ty 

It  (p„u'  + u„p')  * (p„V  ♦ v„p')  + p0  . 0 (9) 
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Equation  of  State 
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(16) 


It  is  noted  that  at  this  point  the  disturbance  field  may  be  non-adiabatic. 

It  was  mentioned  in  the  Introduction  that  the  normal  disturbance 
pressure  variation  is  an  important  mechanism  of  the  lateral  stability  of 
the  flow.  This  rises  from  the  fact  that  the  normal  and  spanwise  momentum 
equations  are  strongly  coupled  through  the  pressure  gradient  terms  even 
though  5p/5y  may  be  small.  However,  at  the  same  time  we  can  neglect  the 
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effect  of  the  pressure  disturbance  in  the  state  equation.  To  see  this, 
consider  the  normal  momentum  equation  estimating  the  orders  of  magnitude 
of  all  entering  terms.  Upon  integration  with  respect  with  y,  the  following 
result  is  obtained  for  the  order  of  magnitude  of  the  pressure: 


- terms  of  order 
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At  the  same  time,  from  the  relation  (16)  we  obtain 
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H'  h-l)Uo  ul 
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Comparing  p'/P0  and  T'/T  (to  simplify,  we  consider  the  adiabatic  case 

H1  = 0 without  loss  of  generality),  it  is  seen  that  when  the  Reynolds 

number  is  large  (6Q/«,  <<  1)  either  p'/P0  <<  T'/T0  or,  if  u'  = 0,  both 

p'/P0  and  T'/To  can  be  neglected  and  hence  p'/p0  - 0-  Both  cases  are 

included  in  the  simplified  state  equation 

r 

M0  0 

We  note  that  this  simplification  cannot  be  used  in  the  region  of  the 

stagnation  point  where  6 /i  >1. 
r oo 

The  pressure  can  be  eliminated  between  the  normal  and  spanwise  momen- 
tum by  cross  derivation  and  subtraction.  As  a result  we  obtain  the  equation 
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We  can  see  now  that  upon  the  assumption  6'  /X  « 1 used  together 
with  the  boundary  layer  approximations,  the  terms  in  the  right  hand  side 
of  the  above  equation  can  be  neglected.  As  a consequence  we  obtain  finally 
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2.3  Assumed  Form  of  Disturbances  and  Resulting  Simplified  Equations 


The  equations  established  above  are  generally  valid  irrespective  of 
the  specific  spanwise  or  streamwise  variation  of  the  disturbances.  Now 
we  specialize  them  to  the  spanwise  periodic  vortex-like  disturbance  field 
mentioned  previously  by  postulating  the  following  forms  for  u',  v',  p',  p', 
T1  near  the  separation  point: 

2 

u'  = u'ty)cos  az  + 0(x-x  ) 


v1  = v^yjcos  az  + 0(x-xs)‘ 

w1  = w^(y)sin  az  + 0(x-xs) 

p‘  = P^(y)cos  az  + 0(x-xs) 

P1  = PgiylCOS  aZ  + 0(x-xs) 

T'  = T^(y)cos  az  + 0(x-xs) 


(19) 
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where  a = 2n/\  is  the  spanwise  wave  number  and  the  subscript  s refers  to 
the  nominal  2-D  separation  point.  These  formulas  include  also  the  require- 
ment that  the  perturbations  undergo  a streamwise  maximum  near  the  separa- 
tion point.  A schematic  illustrating  the  associated  flow  pattern  is  shown 
in  Fig.  1 . 

The  assumed  solution  (15)  reduces  the  number  of  independent  variables 
from  three  to  one,  thus  converting  the  partial  differential  equations  into 
a set  of  ordinary  differential  equations  governing  the  normal  variation  of 
the  disturbances.  At  the  same  time  the  streamwise  gradient  of  the  pres- 
sure disturbance  is  eliminated  from  the  streamwise  momentum  equation. 

Using  the  equations  (14'),  (15),  and  (16)  the  density,  viscosity,  and 
enthalpy  disturbances  can  also  be  eliminated.  Consequently  we  ultimately 

obtain  the  following  set  of  equations  for  u v w and  T 1 : 

s s s s 


3p0  . 
us 


11 


1 


3Ho  u , , . , 3[ho  <V/T0>  * Vs') 

ST  US  * P0v0  3Y 


3H  / 9H  3H 

+ P0  3TVs'  S (Uo  3X2  + P0V0  iT 


8 v T ' 


o \ 

+ c°“°  TF) 


s + u !!!o  s(Ts'/To>+ 
sy-  / To  00  SV  SV 


+ PnVr 


(Ts'/To>  + %us'J 
3V' 


(23) 


where  dY  = p dy  and  h = c T . 

ro  o p o 

We  note  that  when  the  basic  flow  is  adiabatic  (H  = const.)  the 

o 

equation  ( 23)  admits  the  solution  H1  = hQ(Ts 1 /TQ)  + uqus 1 = 0 which 
means  that  the  3-D  disturbed  flow  is  also  adiabatic  (insulating  wall). 


2.4  Disturbance  Equations  in  Terms  of  Levy-Lees  Transformed  Coordinates 


As  will  be  shown  in  the  next  chapter,  Levy-Lees  transformed  coordin- 
ates (including  both  Howarth-Dorodnitsyn  and  Mangier  transformations)  are 
more  convenient  for  the  description  of  the  basic  two-dimensional  flow  near 
the  separation  point  and  hence  for  the  expression  of  the  coefficients 
entering  the  disturbance  equations. 

In  terms  of  Levy-Lees  transformed  coordinates 

. ue  Jo  pody 
25 


= fp  p u 

•'o  w w 


e dx  , 


n = 


(24) 


The  basic  two-dimensional  flow  is  described  by  the  well  known  i ncompressible- 
14- 

like  equations  : 


33f 


32f„  3f  2 Sf  s2f  3f  32f 

rr  * f o-!2  ■ ^0-(^)  ] ♦ f 

o n <?n  on 


3n 


+ f- 


3n 


3f  3g„ 
2f  ( — - — 2. 
* V 3n  35 


3f  3g 
0 a0 

35  3n 


(26) 
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where  3fQ/3n  = uQ/ue,  gQ  - HQ/He  and 


A = - 2 


d In  u H 
e_  _e 

d In  c he 


2;  He 

o „ u 2 dx  ^ 

*\rw  e 


(27) 


Using  the  new  coordinates  and  notations,  the  basic  flow  parameters 
can  be  expressed  by  the  following  relations 

df 


where 


p y u r 

j v - - 2C 

00  L 


T = T * 
o e 


r nJ__  I <)zl 

l>  ' 2c  T * \dn  / J 


37T. 


P e 

, . A . i h«  d 1n  ue 
o ' He  - 'c  H„  d In  { 


= J * 

e 


d In  £ 


(28) 

(29) 

(30) 


X = °e = idL  m 2 ^ 

X 2cpTe*  2 e ^ 


(31) 


Note  that  the  flow  outside  the  boundary  layer  is  isentropic  so  we  have 


H = H , 

6 00 


Te 

T * 
e 


T * = T * 

6 oo 

T / PQ  \ 

oo  i e 


T *\P 

<30  \ OO 


T .*(Cpe)  Y 


III 


duQ  , dC 

„ §L  = _ RT  J 

ue  dx  e C dx 

pe 


where  Cpg  = Pe/PM  and  R is  the  gas  constant. 


(32) 

(33) 

(34) 


We  assume  at  this  point  that  the  functions  fQ  and  gQ  are  regular  at 
the  separation  location*  so  we  can  write 


* This  problem  will  be  discussed  in  more  detail  in  the  next  chapter. 
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2.5  Boundary  Conditions  for  Disturbances  and  Resulting  Eigenvalue  Problem 


The  disturbance  equations  must  be  accompanied  by  the  appropriate 
boundary  conditions  at  the  wall  and  far  from  the  wall.  The  physical  con- 
ditions of  no-slip  and  zero  mass  transfer  (blowing  or  suction)  require 
that  U$ ’ (0)  = Vs ' (0)  = W$'(0)  =0  at  the  wall.  As  for  the  temperature, 
the  problem  is  much  more  complicated  by  the  "spanwise  thermal  response" 

(see  for  example  [ 7])  of  the  solid  body  to  the  thermal  disturbances  in  the 
boundary  layer.  In  the  present  work  we  shall  assume  that  the  surface 
remains  at  the  same  given  temperature  t . Consequently,  we  put  T '(0)  = 0. 

W b 

On  the  other  hand,  far  from  the  wall  and  well  outside  the  boundary  layer 
we  assume  that  the  disturbances  vanish  asymptotically  and  we  put 

us'(«)  = vs'(»)  = ws'(»)  »ts'(«)  =0 

In  practice,  these  conditions  are  implemented  by  constructing  inviscid 
asymptotic  solutions  (far  outside  the  boundary  layer)  which  satisfy  them, 
and  then  imposing  a numerical  matching  of  the  inner  solution  with  the 
asymptotic  values  at  some  finite  but  large  value  of  n. 

The  set  of  differential  equations  (39)  - (42)  along  with  the  afore- 
mentioned boundary  conditions,  when  supplemented  by  the  solution  for  the 
basic  2-D  separating  boundary  layer  profile  functions  and  associated  pres- 
sure gradient  parameters  for  a given  wall  temperature,  constitute  a well- 
posed  (if  rather  complicated)  eigenvalue  problem  for  the  assumed  spanwise 
disturbance  vortices.  Since  the  basic  flow  possesses  a definite  physical 
scale  (namely  the  boundary  layer  thickness  6Q),we  might  expect  to  find  a 
discrete  set  of  eigenvalues  a with  a dominant  fundamental  wavelength  X ~ 
0(6q).  This  similar  to  Gortler's  problem  but  not  the  reattachment  distur- 
bance problem  where  a continuous  spectrum  of  eigenvalues  was  obtained  in 
an  infinitely-wide  incoming  flow. 
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3.  BASIC  UNDISTURBED  2-D  FLOW  DESCRIPTION 
NEAR  A SEPARATION  POINT  (ADIABATIC  CASE) 


3. 1 General  Considerations 

The  problem  of  separation  of  the  laminar  boundary  layer  has  been 
considered  from  the  beginning  as  one  of  the  most  important  for  the  pre- 
diction of  viscous  flows.  The  main  difficulty  arises  from  the  fact  that 
separation  is  the  result  of  the  interaction  between  the  viscous  flow  in 
the  boundary  layer  and  the  external  inviscid  flow.  It  is  well  known  that 
the  boundary  layer  calculation  requires  the  pressure  variation  along  the 
edge,  as  given  by  the  external  inviscid  flow.  In  the  same  time,  since  the 
displacement  thickness  effect  modifies  the  effective  body  shape,  the 
boundary  layer  influences  the  inviscid  flow  and  thus  changes  the  pressure 
distribution.  This  coupled  interaction  becomes  very  strong  near  the  separa- 
tion point. 

In  the  classical  approach,  the  pressure  variation  at  the  edge  is  given 

as  a polynomial  function,  disregarding  this  interactive  displacement  effect 

of  the  boundary  layer.  In  this  case,  as  shown  by  many  authors,  the  boundary 

layer  equation  solution  exhibits  a strong  singularity  in  the  streamwise 

behaviour  of  the  wall  friction  t = (du/dy)  at  the  separation  point.  For 

example,  in  the  case  of  the  flow  past  a plane  with  a given  adverse  pressure 

10 

gradient  of  the  form  dp/dx  = 1 + (x$  - x),  Goldstein  finds  that  du/dy 

1 /2 

behaves  like  (xs  - x)  (Fig.  7a)  and  the  flow  upstream  cannot  be  linked 
up  with  the  flow  downstream  separation  point. 

As  also  shown  by  Goldstein,  it  is  possible  for  a certain  prescribed 
pressure  distribution  that  the  singularity  not  appear  and  (du/dy)wa^  = 0 
for  all  x (e.g.,  the  Falkner-Skan  case).  Goldstein  further  posed  the 
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question  whether  there  might  be  other  special  pressure  distributions  for 
which  the  singularity  does  not  appear.  It  seems  that  Meksyn  gives  an 
affirmative  answer:  using  a different  method  of  analysis,  he  shows  that 

the  singularity  does  not  occur  when  approaching  the  separation  point,  if 
the  streanwise  pressure  variation  corresponds  to  the  real  interactive  one, 
as  experimentally  determined.  The  behaviour  of  (du/dy)  in  such  a case 

Wd  I I 
o 

is  illustrated  in  Fig.  7b  ( as  determined  by  Sfeir  for  the  case  of  the 
flow  in  a corner),  and  we  can  see  that  such  a variation,  although  diffi- 
cult to  be  represented  by  series  expansions,  does  not  exhibit  any  singu- 
larity at  the  separation  point.  The  fact  that  Meksyn  obtains  good  results 
proves  that  the  flow  is  very  sensitive  to  the  pressure  distribution  which 
can  be  considered  reponsible  for  the  occurrance  of  the  singularity,  when 
arbitrarily  given  instead  of  being  allowed  to  float  free  as  part  of  a 
viscous-inviscid  interactive  solution. 

As  also  reported  by  Meksyn,  he  faced  some  difficulties  in  the  compu- 
tations not  in  the  separation  point  but  In  an  upstream  station  close  to 
the  separation  point,  perhaps  in  the  region  of  the  rapid  variation  of 
( du/dy ) wa 1 1 (Fig.  7b). 

In  the  following  Meksyn 's  method  will  be  discussed  in  order  to  adapt 
it  for  the  calculation  of  the  basic  flow  parameters  entering  the  distur- 
bance equations. 


3.2  Adaption  of  Meksyn 's  Method  for  Two-Dimensional  Separation 


If  the  wall  is  insulating  the  energy  equation  ,(26)  admits  the  solu- 
tion gQ  h 1 corresponding  to  the  adiabatic  flow  (Hq  = He  = A = const.). 

The  dynamic  equation  (25)  simplifies,  becoming 


+ f 


A 
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ft,  A\ 

8C  s„2  / 


(45) 


which  Cuincides  in  form  with  the  equation  of  invompressible  laminar 
boundary  layer  studied  by  Meksyn^ . A short  outline  of  his  method  of 
solutioning  in  given  in  Appendix  B. 

In  the  following  we  will  direct  our  attention  to  some  aspects  related 
to  the  use  of  Meksyn's  method  for  the  calculation  of  the  coefficients 
entering  the  disturbance  equations.  As  we  can  see  by  examining  the  equa- 
tions (39),  (40),  (41),  and  (42),  these  coefficients  contain  fQ,  af^/a^, 
3f0/3n.  32f0/3£3n,  32fQ/3n2. 

Using  the  formula  (B-4),  we  can  write  the  following  relations 

no  n d 

f = f dn  /e  “F  <j>dn  = n/  e ”F  <j>  dn  - fe  -Fn4>dn  = nl  - I, 
oo  o o oi 


sT  = nfe  ’F  (ft  ‘ h *)dn  ‘ / 6 _F  n(lt  ‘ n tydr]~ 


T\lr  - i 


0£ 


32f 


H3n  * (e‘F(l{  -I?*)""  ‘ 1 


C£ 


u 

(46) 


where 


2 

= n /e'F  (<t>  + n dn  = 7T 


af  J _F 

I0  = -J  e ^ 4>  dn  (see  also  B-7) 


'°5  ’ { e’F  (ft  ' H *)  d" 


1,  ■ i-f 
0 


1 - j e n dn 


(47) 
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r 


lu 


I3  = / e F (<f>  + n§^-  - n<K)  dn 
0 

and  fQ,  F,  41  are  given  by  (B-3)  and  (B-5). 

2 2 2 

As  we  can  see,  the  calculation  of  fQ,  9fQ/9ri,  9fQ/9£,  9 fQ/9£;9n,  9f  /9n 
reduces  to  the  calculation  of  integrals  IQ,  IQ^,  I-j,  1^  , I3>  all  being  of 
the  form 


/ e "F(c’n)  *(5.n)  dn 


(48) 


where  \p  is  a slowly  varying  function  which  can  be  expanded  in  powers  of 
n (the  coefficients  being  functions  of  c) 


* 


= £ 


C n 


n=0 


(49) 


To  solve  the  integral  (48)  we  use  the  method  of  "steepest  descent"  due 
to  Debye16  and  widely  applied  by  Meksyn  to  calculate  similar  integrals  (see 
also  Appendix  B).  This  method  is  based  on  the  idea  that  the  value  of  the 
integral  (48)  comes  mainly  from  the  vicinity  of  the  stationary  point  of  F, 
which  in  our  case  is  n = 0 (because  9F/9n  = 0 at  n = 0). 


r( 


The  integral  (48)  is  expressed  in  terms  of  incomplete  Gamma  functions 
n+1 


:) 


/ e T ip  dn  = lim  £ d r(  - , x ) 


n+1 


e-*l  n=0 


(50) 


where 


t = F(c,n)  = £ 


n=4 


n-1  n-4 

~nT  n 


(51) 


and  dn's  are  functions  of  a^  (see  Appendix  B)  and  cn>  being  given  in 
Appendix  C. 

To  improve  the  convergence  of  the  series  (50),  we  may  use  Euler's 
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— L I 


Ml 


transformation  which  yields 


J e‘T  ♦ d„  = E P„  2-(n+1  > 
o n=0 


nllC  I C 

P0  = do  r(1/4»  T)»  • • • » Pn  = d0  r(1/4>  +(l)  dl  r^1/2’  T) 


+ f 2 ) d2  r^3/4,  T)  + • • • + dn  r(  — , t) 


/n\  _ n (n- 

w 


[n-m+1 ' 


We  note  that  the  expression  (51)  relating  t and  n can  be  used  only  for 

n < 1.  But  in  the  same  time,  we  can  use  the  following  relation 

n n n 2 

t = f dn  f dn  f e if  dn  = S—  I - n^l  + 2 *2  (53) 


ooo 


where  1 2 = / e"F  n2  <)>  dn  is  also  of  the  form  (48). 
L o 


Because  IQ,  1^  I2  are  calculated  in  terms  of  t,  the  equation  (53) 


gives  n in  respect  with  t 


II  M1!  + 2I0(l/2  I2-t) 


To  find  t by  using  (54)  when  n is  given,  an  iterative  procedure  must  be 
used. 


3.3  Procedure  to  Find  the  Separation  Point  Location  g<;  and  the  Values  of 


(da/dc)c  and  (d2a/dr)c 


As  shown  in  Appendix  B,  fQU,n)  is  expressed  in  terms  of  £,  A(?)»  and 
aU),  where  aU)  is  a given  function  representing  the  pressure  streamwise 
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gradient  while  a(0  must  be  determined  from  the  condition  3fQ/9n  = 1 
at  n = 00  [see  (B-71  )J. 

The  location  of  the  separation  point  can  be  determined  from  the 
condition  (A-l'"),  i.e.,  a(?s)  = 0.  In  this  case  the  relation  (B-71) 


gives 


lim  0 • e + 1/4(24A)1/2  r(l/2)  e2  + 0 

F+l  L 


3 4 5 

+ 0 . e + 0 • e3 


(2 iV/2  (4~15A)A  + 10g  i r,. 

\A  ) 56 

24l_/da\r(2)  8 | 

25  A \d5/  1 ‘ * • • J 


T (3/2 ) e + 0 • 


As  we  can  see  the  first  seven  terms  do  not  contain  da/d^,  this 
appearing  only  in  the  8th  term. 

On  the  other  hand,  by  derivating  the  relation  (B-71)  in  respect  with 


S we  obtain  at  the  separation  location 


1/24 V/4  da 

4\A  / d£ 


r(l/4)  e + 


W) 


^ r(l/2)  e2  + 0 • e3 


* l(l)2  r<’>  * (in 


1/4 r (-141  A2  - 66  A + 94c  ) 


d c ' da 
d£ 


r ( 5/4 ) 


5i  d_r/24\1/2  « - isa)a  ♦ m 

J + d?Lw 


lr  ( 3/2 ) 


- -1  ♦ [&§ 

■ 139^  at  > ♦ ]r(2)  *8  ♦•■•}*  0 (56) 

2 2 

We  can  see  that  the  first  seven  terms  do  not  contain  d a/d£  which  appears 
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only  in  the  eight  term. 


2 2 

To  find  ?s,  (da/d?)s>  and  (da/d?  ) we  can  use  the  following 


procedure: 

1.  Using  the  first  seven  terms  in  Eq.  (55)  we  find  the  separation 
point  ?c 


2.  Using  ? determined  above  and  considering  only  the  first  seven 


terms  in  (56),  we  find  (da/d?), 


3.  Using  (da/d?)  in  the  eight  term  of  (55)  we  recalculate  ? 

^ b 


2 2 

4.  Using  (da/d?)  and  the  new  value  for  ? we  determine  (d  a/d?  ) 

J b b 


from  (56). 

As  a general  remark,  to  solve  the  equations  (55)  and  (56)  for  the 
above  situations  an  iterative  procedure  must  be  used,  starting  from  a point 
? = ?-|  and  approaching  the  separation  point  ? = ?s  where  the  equations  (55) 
and  (56)  are  satisfied.  In  all  cases  the  Euler's  transformation  must  be 
used  to  sum  the  series. 
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4 CONCLUDING  DISCUSSION 


The  theoretical  developments  outlined  in  this  report  represent  the 
background  for  the  investigation  of  the  spanwise  disturbed  two-dimensional 
flow  near  a separation  location.  The  quasi-singular  behaviour  of  the 
basic  flow  near  the  separation,  requiring  a correct,  experimentally  deter- 
mined pressure  distribution  at  the  edge  of  the  boundary  layer,  represents 
the  main  difficulty  of  the  entire  problem.  The  further  steps  towards  a 
final  solution  of  the  problem  are 

1.  Numerical  computation  of  the  basic  flow  parameters  using  the 
formulas  outlined  in  Chapter  3.  In  this  respect,  the  experimental  results 

g 

of  Sfeir  for  the  two-dimensional  separated  flow  in  a corner  can  be 
chosen  as  a concrete  situation.  Due  to  the  simplicity  of  pressure  varia- 
tion for  this  case,  we  can  try  to  approximate  it  by  analytical  formulas  of 
the  form 

p„  = p for  x < x, 

= 1 + A [1  - e‘a(x_xl)]n  for  < x < x2 

' oo 

where  x-j  and  delimit  the  range  of  rapid  pressure  variation  while  A,  a, 
and  n are  parameters  which  can  be  adjusted  in  order  to  obtain  a good  approxi- 
mation for  the  experimental  curve. 

2.  Analytical  and  numerical  study  of  the  eigenvalue  problem  in  order 
to  establish  the  range  of  a for  which  the  disturbance  equations  admit  non- 
trivial solutions. Taking  into  account  the  complexity  of  the  disturbance 
equations  a numerical  computation  seems  to  be  the  only  possibility.  The 
solution  for  the  basic  flow  outlined  in  Chapter  3 suits  only  for  the 


r " — ^ 

adiabatic  case.  An  analytical  treatment  can  be  tried  only  for  an 
incompressible  case. 
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APPENDIX  A 


APPROXIMATE  EXPRESSIONS  FOR  BOUNDARY  LAYER 
INCLINATION  & CURVATURE  AT  SEPARATION 

We  consider  a steady  two-dimensional  high  Reynolds  number  boundary 
layer  flow  without  surface  mass  transfer  and  examine  by  appropriate 
limiting  processes  the  behavior  near  the  surface  at  the  separation  point 
x = x$  C(9u/3y)w>s  = 0]. 

Consider  first  the  streamline  inclination  of  the  flow.  Approaching 
the  surface  at  any  arbitrary  x < x$,  a direct  calculation  of  the  stream- 
line slope  lim  (v/u)  is  indeterminate  since  both  u and  v vanish  at  y = 0; 
y+0 

however,  application  of  L'Hopital's  rule  plus  the  use  of  the  continuity 
equation  resolves  this  difficulty,  giving 

SS-S* 


= lim 
y-+0 


_ -1  9(,pu)  _ 


9X 


3u/9y 


9p 

JLY_ 


(A-l) 


which  always  vanishes  ( as  it  physically  should)  at  the  surface  u = v = 

9u/9x  = 0 for  all  x < x where  (9u/9y)  t 0.  Approaching  the  separation 

s w 

point  where  (3u/9y),,  -*■  0,  however,  the  right  side  of  A-l  becomes  indeter- 
w 

minate  so  we  must  again  apply  L'Hopital  to  obtain 
lim  (v/u) 
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(A-2) 


5 • [»x  (ay)/ sy  (w)]„al1,x 

which  holds  for  an  arbitrary  streamwise  pressure  gradient,  degree  of  com- 
pressibility and  turbulent  eddy  viscosity  model.  But  for  any  such  flow 
we  can  write  t = (u  + Uy)  (3u/3y)  while  approaching  the  wall  3i/3y  = 

3p  /3x ; using  these  in  (A-2)  and  assuming  that  any  turbulence  damps  out  at 
the  wall  (uyw  = 0),  we  obtain 


1 im 
y-K) 
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•o  W, 


" (3rw/3x)s 

(V9X)S 


(A-3) 


a result  again  valid  for  a perfectly  general  flow.  Equation  A-3  defines 
the  slope  of  the  separation  streamline  leaving  the  surface  (Fig.  4). 

Now  consider  the  average  curvature  of  the  flow  passing  above  this 
separation  streamline.  Assuming  it  is  not  large,  we  can  write  the  associ- 
ated radius  of  curvature  to  good  approximation  as 

rAV  - “ax  ■ u + u u (A_4) 

Approaching  the  surface  (which  is  presumed  flat)  upstream  of  the 
separation  point,  this  expression  is  indeterminate;  then  once  more  using 
L'Hopital  gives 
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at  x = x. 


Resolving  the  latter  indeterminate  form,  we  get 
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= lim  3x‘ 
y+0 
x->x. 


*)•©>(§) 


3/  3u> 

3y  1 3y, 


which  upon  differentiations  of  the  formulae  T = (u  + Uj)  (3u/3y)  and 
3x/3y  = 3p/3x  and  assuming  (3ux/3y)  = 0 ultimately  yields  the  following 

i W j s 


expression: 


(f\\  M ?>1  3_KY^ 
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2 2 

Since  in  most  applications  3 p /3x  and  3u  /3x  * 3T./3x  are  very  small 

w w w 

if  not  in  fact  zero  at  x$  while  the  corresponding  heat  transfer  and  hence 
(3u/3y)s  is  also  usually  small,  it  is  sufficient  to  neglect  all  but  the 


first  term  in  the  numerator  and  get 


(3PW/3X)S 


(3  V8X  >s 


(A-6) 


We  observe  from  (A-6)  that  Rftv  $ is  directly  proportional  to  the 
curvature  of  the  t (x)  plot  and  that  it  also  decreases  as  the  separation 
angle  increases  (i.e.,  in  agreement  with  physical  expectation,  the  curva- 
ture vanishes  as  this  angle  -+  0). 


APPENDIX  B 


MEKSYN'S  SOLUTION  FOR  THE  TWO-DIMENSIONAL 
FLOW  IN  THE  BOUNDARY  LAYER  NEAR 


SEPARATION 


The  problem  is  to  find  the  solution  of  the  equation  (39),  the 
boundary  conditions  being 


= 0 at  n = 0 


(B-T) 


= 1 at  n = 


(B-l ' ' ) 


At  the  separation  point  there  is  an  additional  condition 


= 0 at  n = 0 


(B-T") 


As  shown  by  Meksyn,  in  the  region  close  to  the  wall  (n  < 1),  f (?,n) 
can  be  expanded  in  powers  of  n,  the  coefficients  being  functions  of  ?, 
A(£),  and  a ( C ) = (32fQ/3n2)n=0 
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In  the  neighbourhood  of  the  separation  point  Meksyn  gives  the  following 

expressions  for  a coefficients: 
n 

a2  = 0 [hence  fQ(?,n)  starts  with  n3,  see  (B-l"')]; 
a3  = A;  a4  = 0;  ag  = 2?  a (da/d?); 

ag  = - (4  + 6A)Aa  + 4?a  (dA/d?);  (B_3) 

a-,  = - (4  + 6a)a2  + 4?A  (dA/d?); 

ag  = - 4?2a  (da/d?)2; 

ag  = - 2(10  + 42A)?Aa  0 + 2?2[-26a  |}A  4|  + 18  ^ + 18Aa^-|] 


31 


The  solution  of  the  equation  (B-l)  can  also  be  expressed  in  the 


form 


32f 


= /fn1  n 


(B-4) 
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where  <f>( C » n)  is  a function  whose  expansion  in  powers  of  n near  the  wall  is 

given  by 


4>U»n) 
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(B-5) 


In  the  region  of  separation  bn 1 s have  the  following  expressions 


da 


bQ  = a;  b]  = A;  b 2 =0;  b3  = 2^a  ; 


dA 


b4  = -3(1  + 2a) Aa  + 4c  ^ a ; 


b5  = (1  - 6A)A2  + 4CA  ^ ; 


(B-6) 


b6  = -45  a 


Meksyn  demonstrates  that  the  relation  (B-4)  also  describes  the  asymp- 


p p 

totic  behaviour  of  (3  f /3n  ) far  from  the  wall,  which  shows  that  <j>(S»n) 


must  be  considered  a slowly  varying  function. 

Although  the  expansion  (B-5)  of  4>  is  divergent  for  n > 1,  the  relation 
(A-4)  can  still  be  used  to  approximate 

n 


af 


= / e F <|>(C  ,n)  dn 


because  F has  n = 0 as  a stationary  point,  and  as  a consequence,  the  value 
of  the  integral  comes  mainly  from  the  points  near  the  wall  where  F changes 
very  rapidly. 


Using  this  idea,  Meksyn  gives  the  following  expression  for  9fQ/3n  in 


-a- 
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n+ 1 

terms  of  t = F(n)  and  incomplete  Gamma  functions  r(  pp  * t) 
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&fo  f -F  . . v a n+1  \ n+l 

inr  ■ J e * d"  ' ’>?  dn  r(  T ' t!  c 

o g-m  n=0 

where  dn's  have  the  following  expressions  near  the  separation  point 

do  ■?(fi),/4  a ; d,  ' i (24»>'/2-  d2  = 0: 
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The  unknown  function  a(^)  entering  the  above  relations  is  found  by  using 

the  boundary  condition  (B-l")  which  furnishes  the  following  differential 


equation 


lim  £ dn  r(  ) e"+1  - 1 
e-*l  n=0  n 4 


(B-71) 


where  r(  pp  ) are  complete  Gamma  functions. 

To  improve  the  convergence  of  the  above  series  when  e-*l , use  must  be 
made  of  the  well-known  Euler's  transformation  [see  also  (52)],  this  giving 


lim  L p e n+1  = 1 
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• t 


r(l/4),  . . 


pn  = dQ  r(i/4)  + (Jj  d1  r(l/2)  + d2  r(3/4) 


Usually,  the  first  six  terms  in  ( B-7 1 ' ) are  sufficient. 


APPENDIX  C 


d COEFFICIENTS  ENTERING  THE  RELATION  (50) 
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^ REGION  "A" 


Reattachment  Flow  Configuration  (Schematic) 


Spanwise-Periodic  Disturbances  in  Reattaching  Flows  linger) 
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Fig.  3 Detailed  Surface  Flow  Pattern  Near  Ramp-Induced  Separation  Line 

(Settles,  et  al) 
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Fig.  4 Schematic  of  Assumed  3-D  Disturbance  Field  in  Separating  Flow  Region 


Schematic  Illustration  of  Separation  Region  Inclination  and  Curvature  Estimation 


Pressure  Distribution  Along  a Ramp-Induced  Separation  Region  (Steir) 
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